Periodic Anderson model with d-f interaction 
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We investigate an extended version of the periodic Anderson model where an interaction is 
switched on between the doubly occupied d- and /-sites. We perform variational calculations us- 
ing the Gutzwiller trial wave function. We calculate the /-level occupancy as a function of the 
/-level energy with different interaction strengths, ft is shown that the region of valence transition 
is sharpened due to the new interaction. 
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I. INTRODUCTION 

The heavy-fermion phenomenon, which is observable 
in many rare-earth materials is still an active field of re- 
searched area. The simplest model, which can account 
for these phenomena is the periodic Anderson model. For 
a review on this topic see [lH3j. The model in its sim- 
plest form describes a hybridization (v) between a wide 
conduction band and correlated /-electrons localized at 
lattice sites. Despite its simple form there is no gen- 
eral solution to this problem, exact solutions have been 
obtained only in some special cases Nonperturba- 
tive approaches like variational methods [5-9] are appli- 
cable even if the Coulomb interaction (£7/) between the 
/-electrons is large, though we have to use some uncon- 
trolled approximation when calculating the expectation 
values. In a recent paper 10] it was shown that the 
repulsion between the d- and /-electrons plays an im- 
portant role in valence transition. It has been argued 
that the sudden change of valence of Ce can account 
for the presence of a sharp peak in the transition tem- 
perature of some superconducting Ce based compounds. 
The effects of this new interaction have been examined 
with dynamical mean field theory HM3, slave-boson 



method [ljj [H| i variational Monte Carlo technique [15 1 
and projector-based renormalization approach [Ig. It 
has been found in these papers that intermediate valent 
region narrows in the presence of d-f coupling. 

In this paper we examine the ground state properties of 
this extended model using the Gutzwiller method. How- 
ever, the conventional d-f interaction 
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is too difficult to handle within the framework of this 
method, so we consider a modified d-f interaction 
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where electrons repel each other only when site j is fully 
occupied. The reason for this modification is that it can 
be treated easier, since much less electron configurations 
need to be taken into account. 



II. VARIATIONAL CALCULATION 

We summarize the main steps of the performed vari- 
ational calculation following Ref. [13] ■ For simplicity 
in the present work we restrict ourselves to a nonmag- 
netic half-filled band, although the variational procedure 
can be carried out without these assumptions. Assuming 
that there are iV-f = N± up- and down-spin electrons in 
an arbitrary dimensional lattice with L lattice sites, we 
consider the following Hamiltonian: 

h = E^( fe )4,A- - «E i&A* + 
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where W j a = /j )<7 /j>, and h d j a = dj a dj tfT . The symbol 
k denotes the wave vector, and j the site index. The 
width of the d-band is W . The trial wave function of 
Gutzwiller-type is expressed as 

i*) =^nn^/L+^ijio), (4) 
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where the mixing amplitudes Uk/vk are also variational 
parameters, and Pq and Pq are the Gutzwiller projec- 
tors, which are written down as follows: 

Pi = W - (! - n^n^h^}, (5) 
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Pg = Illl-il-VfWn*}^ (6) 
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where the variational parameters r/df and rjf are con- 
trolled by Udf and Uf, respectively. We use the so- 
called Gutzwiller approximation to evaluate the expec- 
tation values of the Hamiltonian. Optimizing the ground 
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state energy density with respect to the mixing ampli- 
tudes, we obtain the following expression for that: 

fceFS L 

+0/ -e/W +UdfVdf + U f Uf, (7) 
where is the average number of /-electrons per site, 



i.e. 



71/ = 



(S) 



v j denotes the density of the doubly occupied /-sites and 
Vdf is the density of such sites when both the /- and d- 



sites are doubly occupied, e/ is the quasiparticle energy 
level of the /-electron (the effective /-level, its precise 
form is not presented here due to its lengthy form) , v de- 
notes the renormalized hybridized amplitude, v=^JqdqfV, 
and the summation is carried out over all wave numbers, 
fc, of electrons in the Fermi sea (FS). The qd and qf 
are the kinetic energy renormalization factors for the d- 
electrons and /-electrons respectively and are functions 
of n (the total numbers of electrons per site), nf,v d (the 
density of the doubly occupied d-sites), Vf, and Vdf- The 
determination of these quantities is the main task of the 
Gutzwiller method. Their forms in the present model are 
much more complicated than in those models which do 
not contain the interaction U d f. The results are written 
in the following complete square (Gutzwiller-likc) form: 
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where n is the total number of electrons per site, i.e. 

n = T (Si,a + Ej,a ^,<r) ■ lt is wortn mentioning 

that without Udf, qf depends only on nf and Vf p7| . 
Since there is no direct hopping between the different 
/-sites, qf appears only in the renormalized hybridized 
amplitude v~^jq d qfv. It should be emphasized that q d 
is not unity, i.e. the width of the <i-band is narrowed, 
though there is no Coulomb repulsion between the d- 
electrons. It is also remarkable that the renormalization 
amplitude of the hybridization still can be written as the 
square root of the product of qf and q d in spite of this 
new interaction. These are the main results of this paper. 



III. RESULTS 

For further calculations we assume that the density of 
states of the <i-band is constant, p(e) — 1/W, when e is 
in the interval [-W/2, W/2]. Using this assumption the 
summation over the wave vectors in Eq. ([7]) can still be 
performed analytically. The ground state energy density 
thus obtained should be optimized with respect to v d , 
Vf, v d f and nf. After this procedure we arrive at a non- 



linear system of equations for these unknown quantities. 
We solve the obtained equations numerically. Typical 
results are shown in Fig. [TJ One can see that switching 




FIG. 1: The /-level occupancy as a function of the energy e/ 
of the /-level at half-filling, v/W = 0.1875 and U f /W = 5. 
The continuous and dotted curves belong to Udf /W = and 
10, respectively. The inset shows the enlarged view of the 
framed region. 
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on Udf also results in the narrowing of the intermediate 
valent regime of 2> nj >1 as the conventional d-f in- 
teraction (see Eq.((T|)) while no narrowing occurs 
in the region of 1> rif >0, since the d-f interaction we 
used is effective only when there are many doubly oc- 
cupied /- and d-sites. The conventional d-f interaction 
significantly narrows the intermediate valent regime of 
1> 71/ >0, too. 

It is interesting to investigate the q/ renormalization 
factor, since it is related to the effective mass (m* ~ qj ). 
The results are shown in Fig. [5J The same effect can 
be seen here too, furthermore the heavy-fermion regime 
extends to smaller /-level energies. 

In Fig. [3] the renormalization factor for the d-band is 
plotted. As it was outlined in Eq. (fTU|) the d-electrons 
become correlated in the presence of Udf, despite the fact 
that there is no direct d-d Coulomb interaction present 
in the model. 




FIG. 2: The /-level renormalization factor as a function of 
the energy e/ of the /-level at half-filling, v/W — 0.1875 and 
Uf/W — 5. The continuous and dotted curves belong to 
Udf /W — and 10, respectively. 
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FIG. 3: The d-band renormalization factor as a function of 
the energy e/ of the /-level at half-filling, v/W = 0.1875 and 
Uf/W — 5. The continuous and dotted curves belong to 
Udf /W — and 10, respectively. 



IV. CONCLUSIONS 



We have investigated an extended version of the usual 
periodic Anderson model, where an additional repulsive 
interaction among the d and / electrons was taken into 
account. We have shown analytically that the conduc- 
tion band is narrowed owing to the new interaction. We 
have demonstrated that a finite Udf expands the regime 
of heavy-fermion character to a lower value of e/ and nar- 
rows the intermediate valent regime at the same time. It 
is expected that the full d-f interaction (Eq. (JTJ) has an 
even sharper effect. The calculation is in progress. 
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